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ABSTRACT 

We present the results from a detailed ASCA/ROSAT X-ray study of three distant, 
massive cooling flows; Zwicky 3146 (z = 0.291), Abell 1835 {z = 0.252) and E1455+223 
(Zwicky 7160; z = 0.258). Using multiphase models fitted to the ASCA spectra, we 
determine values for the temperature, metallicity, luminosity and cooling rates in 
the clusters. These results are combined with deprojection analyses of the ROSAT 
images to provide detailed constraints on the mass distributions in the systems, and 
on the properties of their cooling flows. The spectral and imaging data identify these 
clusters as the three most-massive cooling flows known, with mass deposition rates of 
~ 1400, 2300, and 1500 Mq yr~^ respectively. We highlight the need for multiphase 
models to consistently model the spectral and imaging X-ray data, and discuss the 
importance of using these models in X-ray determinations of the cluster masses. We 
also present results from an extensive optical study of the clusters and report the 
discovery of a gravitational arc in Abell 1835. The lensing data provide an independent 
constraint on the distribution of mass in the cluster, in good agreement with the 
results from the multiphase (although not single-phase) X-ray analysis. We present 
measurements of the galaxy distributions in Abell 1835 and E1455-I-223 and relate 
these to the distributions of the total mass. The ASCA spectra place firm constraints 
on the column density of intrinsic X-ray absorbing material in the clusters. Abell 1835 
and E1455+223 exhibit large intrinsic column densities (TVh ~ 3 — 4 x 10^^ atom cm~^) 
associated with their cooling flows. These clusters also exhibit significant amounts of 
reddening in the optical spectra of their central galaxies. The data for Zwicky 3146 
indicate lower lower levels of X-ray absorption and reddening. All three clusters exhibit 
excellent alignment between the position angles of their X-ray emission and the optical 
emission from their dominant cluster galaxies. 

Key words: galaxies: clusters: general - cooling flows - intergalactic medium - dark 
matter - gravitational lensing - X-rays: galaxies 



1 INTRODUCTION 

X-ray observations of clusters of galaxies show that in the 
central regions of most clusters the cooling time of the Intr- 
aCluster Medium (ICM) is significantly less than the Hub- 
ble time (e.g. Edge, Stewart & Fabian 1992). The observed 
coofing, which takes place primarily through the emission of 
X-rays, leads to a slow net inflow of material towards the 
cluster centre; a process known as a cooling flow (Fabian 
1994). 

Within the idealized model of an homogeneous cooling 
flow, all of the cooling gas flows to the centre of the cluster 
where it is deposited, having radiated away its thermal en- 



ergy. However, observations of cooling flows show that the 
simple homogeneous model is not correct. Although the X- 
ray surface brightness profiles of clusters with cooling flows 
are substantially more sharply-peaked than those of non 
cooling-flow systems, the emission is not as sharply-peaked 
as it would be in the case of homogeneous flows. Rather 
than all of the cooling gas flowing to the very centres of 
the clusters, the X-ray data show that gas is 'cooling out' 
throughout the central few tens to hundreds of kpc. Typi- 
cally the cooled material is deposited with Af(r) cx r, where 
M{r) is the integrated mass deposition rate within radius r. 
The X-ray data flrmly require that cooling flows are inho- 
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mogeneous with a range of density and temperature phases 
at all radii. 

Spatially resolved X-ray spectroscopy of clusters also 
confirms the the presence of distributed cool (and rapidly 
cooling) gas in cooling flows. The spatial distribution and lu- 
minosity of the cool components determined from the spec- 
tral data are well-matched to values inferred from the X-ray 
images {e.g. Allen et al. 1993, Allen, Fabian & Kneib 1996; 
Fabian et al. 1996; Allen & Fabian 1996). For a detailed 
review of the theory and observations of cooling flows see 
Fabian (1994). 

In this paper we present observations, made with the 
ASCA and ROSAT X-ray astronomy satellites, of three ex- 
ceptionally luminous cluster cooling flows. Two of the sys- 
tems, Zwicky 3146 {z = 0.291) and Abell 1835 {z = 0.252) 
were identified as X-ray luminous clusters during optical- 
foUow up studies to the ROSAT All-Sky Survey (Allen et 
al. 1992a). The third system, E1455+223 (or Zwicky 7160; 
z = 0.258) was identified by Mason et al. (1981) in a follow- 
up to X-ray observations made with Einstein Observatory. 
All three clusters are included in the ROSAT Brightest Clus- 
ter Sample (Ebeling et al. 1996a). 

Combining the high spectral resolution ASCA data with 
ROSAT images we present consistent determinations of the 
temperatures, metallicitics, luminosities and cooling rates 
in the clusters. The data for Zwicky 3146, Abell 1835 and 
E1455-I-223 identify them as the three largest cluster cooling 
flows known to date. We constrain the level of intrinsic X- 
ray absorption in the cooling flows and relate the results to 
measurements of intrinsic reddening in the Central Cluster 
Galaxies (CCGs) of the systems. Results on the distributions 
of X-ray gas, galaxies, and the total gravitating matter in 
the clusters are reported. 

The structure of this paper is £is follows. In Section 2 
we summarize the observations. In Section 3 we present the 
X-ray and optical imaging data and discuss the morpholog- 
ical relationships between the clusters and their CCGs. In 
Section 4 we discuss the spectral analysis of the X-ray data. 
In Section 5 present the results from the deprojection anal- 
yses of the ROSAT images. Section 6 discusses the optical 
properties of the clusters. In Section 7 we discuss some of 
the more important results in detail, and in Section 8 sum- 
marize our conclusions. Throughout this paper, we assume 
Ho=50 kms"^ Mpc"\ = 1 and A = 0. 



2 OBSERVATIONS 

The details of the observations are summarised in Table 1. 
Exposure times for the ASCA data sets arc effective expo- 
sures after standard data screening and cleaning procedures 
have been applied (Day et al. 1995). Hashed numbers fol- 
lowing an instrument name indicate that those observations 
were carried out on more than one date. The ASCA observa- 
tions of Abell 1835 were made in two parts, on consecutive 
days in 1994 July. The ROSAT HRI observations of Zwicky 
3146 were carried out on 2 dates, in 1992 November and 
1993 May, giving a total exposure of 26.0 ks. The HRI ob- 
servations of E1455+223 were carried out in 3 parts, in 1992 
January - 1994 July, with a total exposure time of 14.9 ks. 

The ASCA observation of Zwicky 3146 was carried out 
in 1993 May during the PV stage of the mission. The SIS 



detectors were used in 2-CCD mode with the target posi- 
tioned approximately at the boundary between chips 1 and 
2 in SISO (the nominal pointing position in 2-CCD mode 
during the PV phase). The observations of Abell 1835 and 
E1455+223 were carried out during the AO-1 stage of the 
ASCA program. These clusters were also observed in 2- 
CCD mode, but with the targets positioned more centrally 
in chip 1 of SISO. For a detailed discussion of ASCA ob- 
serving modes and instrument configurations see Day et al. 
(1995). Reduction of the ASCA data was carried out using 
the FTOOLS package. Standard selection and screening cri- 
teria were applied (Day et al. 1995). The ROSAT data were 
analysed using the STARLINK ASTERIX package. 

The optical observations of Zwicky 3146 was carried out 
with the 3.6m telescope at the European Southern Obser- 
vatory (ESO), La Silla, Chile. The ESO Faint Object Spec- 
trograph Camera was used with the TEK 512 x 512 CCD 
(pixel scale 0.61 arcscc). An exposure of 120 sec in the R 
band was made in seeing of ~ 1.5 arcsec. The optical ob- 
servations of Abell 1835 and E1455-I-223 were carried out 
with the 5m Hale Telescope, Palomar, as part of a follow up 
study (Edge et al. 1996) of the most X-ray luminous clusters 
in the ROSAT Brightest Cluster Sample. The COSMIC in- 
strument and TEK 2048 x 2048 chip (pixel scale 0.28 arcsec) 
were used. Exposures of 1000 and 500 s were made in Gunn 
i, 500s and 600s in KC B, and 3000s in KG U, for Abell 1835 
and E1455-I-223 respectively. The seeing was ~ 1.1 arcsec. 
The optical data were reduced and analysed in IRAF. 

Figs. 1-3 show the optical images of the clusters, with 
the ROSAT HRI X-ray contours overlaid. (Details of the 
smoothing algorithms are given in the figure captions.) 



3 MORPHOLOGY ANALYSIS 

The images presented in Figs. 1-3, and the optical and X-ray 
co-ordinates listed in Table 2, demonstrate excellent agree- 
ment between the positions of the CCGs and the positions 
of the peaks of the X-ray emission from the clusters. [Errors 
of < 5 arcscc may be associated with the aspect solutions of 
the HRI data. The astrometry of the optical data is accurate 
to within 1 arcsec] 

The ellipticities and position angles of the X-ray emis- 
sion from the clusters and the optical emission from the 
CCGs have been examined using the ELLIPSE isophote- 
analysis routines in IRAF. The images were re-binned to 
a suitable pixel size (8x8 arcsec^ for the X-ray data and 
0.57x0.57 arcscc'^ for the optical images) and were modelled 
with elliptical isophotes (Jedrzejewski 1987). The elliptici- 
ties, position angles and centroids of the isophotes were free 
parameters in the fits. The results are summarized in Ta- 
ble 3 where we list the mean ellipticities and position angles 
over the range of semi-major axes studied. We find excellent 
agreement between the position angles of the optical (CCG) 
and X-ray (cluster) isophotes. 

The agreement between the position angles of the CCG 
and cluster isophotes, and the coincidence of the CCG co- 
ordinates and the peaks of the cluster X-ray emission, are 
similar to the results from studies of other large cooling- 
fiow clusters at lower redshifts (White et al. 1994; Allen et 
al. 1995; Allen et al. 1996). 
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Table 1. Observation summary 



Cluster 


Instrument 


Observation Date 


Exposure (ks) 


Zwicky 3146 


ASCA SISO 


1993 May 18 


26.3 




ASCA SISl 




31.6 




ASCA GIS2 




30.2 




ASCA GIS3 




30.2 




ROSAT HRI #1 


1992 Nov 27 


15.2 




ROSAT HRI #2 


1993 May 17 


10.8 




ROSAT PSPC 


1993 Nov 13 


8.62 




ESO 3.6m (R) 


1992 Nov 11 


0.12 


Abell 1835 


ASCA SISO #1 


1994 Jul 20 


18.2 




ASCA SISl #1 


55 55 


17.2 




ASCA GIS2 #1 


55 55 


13.0 




ASCA GIS3 #1 




13.0 




ASCA SISO #2 


1994 Jul 21 


8.51 




ASCA SISl #2 




8.22 




ASCA GIS2 #2 




6.55 




ASCA GIS3 #2 




6.56 




ROSAT HRI 


1993 Jan 22 


2.85 




ROSAT PSPC 


1993 Jul 03 


6.18 




Hale 5m (Gunn i) 


1994 Jun 09 


1.00 




Hale 5m (KG B) 


1994 Jun 10 


0.50 




Hale 5m (KG U) 


1994 Jun 09 


3.00 


E1455+223 


ASCA SISO 


1994 Jul 18 


30.5 




ASCA SISl 


55 55 


28.3 




ASCA GiS2 


55 55 


18.3 




ASCA GIS3 




18.3 




ROSAT HRI #1 


1992 Jan 11 


4.09 




ROSAT HRI #2 


1993 Jan 20 


4.23 




ROSAT HRI #3 


1994 Jul 07 


6.57 




Hale 5m (Gunn i) 


1994 Jun 10 


0.50 




Hale 5m (KG B) 


1994 Jun 10 


0.60 




Hale 5m (KG U) 


1994 Jun 10 


3.00 



Notes: Exposure times are effective exijosurcs after all cleaning and correction procedures have been carried out. For Abell 1835 the ASCA 
observations were carried out in 2 parts, yielding total exposures of 26.7, 25.4, 19.6 and 19.6 ks for the SO, SI, G2 and G3 detectors, 
respectively. For Zwicky 3146, two separate ROSAT HRI observations were made, yielding a total exposure of 26.0 ks. For E1455+223 
three HRI observations were carried out providing a total exposure of 14.9 ks. 



Table 2. Target summary 







OPTICAL (J2000.) 


X-RAY (J2000.) 






Cluster 


z 


R.A. 


Dec. 


R.A. 


Dec. 


Fx 


Lx 


Zwicky 3146 


0.2906 


10'^23™39.6'' 


04°11'10" 


10'>23™39.8= 


04°11'11" 


6.6 


2.8 


Abell 1835 


0.2523 


14'>01™02.0" 


02°52'42" 


14hoimoi.9s 


02°52'43" 


12.3 


3.8 


E1455+223 


0.2578 


14h57™15.P 


22°20'31" 


14^57™ 15.0= 


22°20'36" 


3.7 


1.3 



Notes: Redshifts and CCG coordinates (J2000) from Allen et al. (1992). X-ray coordinates denote the position of the X-ray peak 
determined from the HRI data. X-ray fluxes {Fx) in units of 10^^^ ergcm"^ s~^ and luminosities (Lx) in 10*® ergs^^ are determined 
from the ASCA (SO) data. Fluxes are quoted in the 2-10 keV band of the observer. Luminosities are absorption-corrected and are quoted 
in the 2-10 keV rest-frame of the object. 
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Figure 1. The ESO 3.6m R band image of Zwicky 3146 with the ROSAT HRI X-ray contours overlaid. The optical data have a pixel 
scale of 0.61 arcsec and were taken in ~ 1.5 arcscc seeing. The X-ray image has a pixel size of 2 X 2 arcscc^ and has been adaptivcly 
smoothed (Ebeling, White & Rangarajan 1996) to give > 36 count smoothing element"^. Contours are drawn at eight evenly-spaced 
logarithmic intervals between 0.89 and 22.38 ct pixel"^. 



Table 3. Isophote Analysis 



Cluster OPTICAL X-RAY 

pixel range ellipticity RA. pixel range ellipticity P.A. 

Zwicky 3146 0.61 x 0.61 1.2 - 3.7 0.32 ± 0.01 125 ±1 8.0 x 8.0 16 - 24 0.16 ± 0.04 127 ± 8 

Abell 1835 0.57 x 0.57 1.1 -4.5 0.21 ±0.03 145 ±4 8.0 x 8.0 16-56 0.20 ± 0.07 163 ± 16 

E1455-I-223 0.57 x 0.57 1.1 -4.5 0.16 ±0.01 35 ± 2 8.0 x 8.0 16-48 0.20 ± 0.05 40 ± 7 



Notes: A summary of the results from the isophote analysis of the optical (CCG) and X-ray (cluster) data. Columns (2) and (6) give 
the pixel sizes in arcsec of the re-binned optical and X-ray images used in the analyses. Columns (3) and (7) list the range (in arcsec) 
of semi-major axes analysed. Columns (4) and (8) give the mean ellipticities (defined as 1 — b/a where b and a axe the semi-minor and 
semi-major axes respectively) of the optical and X-ray data in the regions analysed. Columns (5) and (9) list the mean position angles 
(PA) in degrees in these same regions. 



4 SPECTRAL ANALYSIS OF THE ASCA DATA 
4.1 Method of Analysis 

SIS spectra were extracted from circular regions centred on 
the positions of the X-ray peaks. The radii of these regions 



were selected to minimize the number of chip boundaries 
crossed (and thereby minimize systematic uncertainties in- 
troduced into the data by such crossings) whilst covering as 
large a region of the clusters as possible. The compromise of 
these considerations lead to the choice of regions for spectral 



© 0000 RAS, MNRAS 000, 000-000 



ASCA and ROSAT observations of distant massive cooling flows 5 



Figure 2. The Hale 5m Gunn i image of Abell 1835 with the ROSAT HRI X-ray contours overlaid. The optical data have a pixel scale of 
0.28 arcscc and were taken in ~ 1.1 arcscc seeing. The X-ray image has a pixel size of 4 X 4 arcsec-^ and has been adaptively smoothed to 
give > 16 count smoothing element"^. Contours are drawn at six evenly-spaced logarithmic intervals between 0.79 and 7.94 ct pixel"^. 



Table 4. Regions included in the spectral analysis 



Cluster 


Detector 


Radius (arcmin/kpc) 


Zwicky 3146 


SISO 


5.0/1620 




SISl 


4.0/1230 




SIS2 


6.0/1940 




SIS3 


6.0/1940 


Abell 1835 


SISO 


4.0/1190 




SISl 


3.0/890 




SIS2 


6.0/1780 




SIS3 


6.0/1780 


E1455H-223 


SISO 


3.8/1140 




SISl 


2.8/840 




SIS2 


6.0/1800 




SIS3 


6.0/1800 
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analysis listed in Table 4. For Abell 1835 and E1455-I-223 the 

spectra were extracted from a single chip in each SIS (chip 1 
in SISO and chip 3 in SISl). For Zwicky 3146, which is cen- 
tred on the boundary between chips 1 and 2 in SISO (chips 
3 and in SISl), the data were extracted across the two 
chips in circular regions bounded by the outer chip edges. 
The GIS spectra used for the analysis of the cluster proper- 
ties were extracted from circular regions of radius 6 arcmin 
(corresponding to ~ 2 Mpc at the redshifts of the clusters) 
again centred on the peak of the X-ray emission from the 
clusters. 

Background subtraction was carried out using the 
'blank sky' observations compiled during the performance 
verification stage of the ASCA mission. (The blank sky ob- 
servations arc compiled from observations of high Galac- 
tic latitude fields free of bright X-ray sources). For X-ray 
sources lying in directions of relatively low Galactic column 
density, like the targets discussed in this paper, the blank 
sky observations provide a reasonable representation of the 
cosmic and instrumental backgrounds in the detectors over 
the energy ranges of interest. 

The modelling of the X-ray spectra has been carried 
out using the XSPEC spectral fitting package (version 8.50; 
Shafer et al. 1991). For the SIS analysis, the 1994 Novem- 
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Figure 3. The Hale 5m Gunn i image of E1455+223 with the ROSAT HRI X-ray contours overlaid. The optical data have a pixel 

scale of 0.28 arcscc and were taken in ~ 1.1 arcscc seeing. The X-ray image has a pixel size of 2 X 2 arcscc^ and has been adaptively 
smoothed to give > 25 count smoothing element"^. Contours are drawn at nine evenly-spaced logarithmic intervals between 0.21 and 
8.32 ct pixel" 1. 



bcr 9 release of the response matrices from GSFC was used. 
Only those counts in pulse height analyser (PHA) channels 
corresponding to energies between 0.5 and 10 keV were 
included in the fits (the energies between which the calibra- 
tion of the SIS is best-understood). For the GIS analysis, 
the 1995 March 6 release of the GSFC response matrices 
was used and only data in the energy range 1 — 10 keV were 
included in the fits. All spectra were grouped before fitting 
to ensure a minimum of 20 counts per PHA channel, thereby 
allowing statistics to be used. 

The X-ray emission from the clusters has been mod- 
elled using the plasma codes of Raymond & Smith (1977; 
with updates incorporated into XSPEC version 8.50) and 
Kaastra & Mewe (1993). The results for the two plasma 
codes show good agreement. For clarity, only the results for 
the Raymond & Smith (hereafter RS) code will be presented 
in detail in this paper although the conclusions drawn may 
equally be applied to the analysis with the Kaastra & Mewe 
code. 

We have modelled the ASCA spectra both by fitting 
the data from the individual detectors independently and by 
combining the data from all 4 detectors. The results from 
the fits to the individual detectors are summarised in Ta^ 
bles 5-7. When combining the data for the different detec- 



tors the temperature, metallicity and column density values 
were linked together. However, the normalizations of both 
the ambient cluster emission and the cooling flow compo- 
nents were allowed to vary independently, due to the differ- 
ent source extraction areas used and residual uncertainties 
in the flux calibration of the instruments. The results from 
the fits to the combined data sets are summarized in Table 
8. 



4.2 The spectral models 

The ASCA spectra were first examined with a simple single- 
phase model consisting of an RS component, to account for 
the X-ray emission from the cluster, and a photoelectric 
absorption component (Morrison & McCammon 1983) nor- 
malized to the equivalent Galactic hydrogen column density 
along the line-of-sight to the cluster. The free parameters 
in this model (hereafter Model A) were the temperature, 
metallicity and emission measure of the X-ray gas. The red- 
shift of the X-ray emission from the cluster was fixed at the 
optically-determined values for the CCGs (Table 2). 

We then examined a second model (Model B) in which 
the absorbing column density was also allowed to be a free 
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Table 5. Spectral Analysis of the ASCA data for Zwicky 3146 



7 





Parameters 


SO 


SI 


G2 


G3 


MODEL A 


kT 
Z 

Nh 
X /uut 


5 6+0-3 
22+° °'' 
0.34 

z4o. 1 1 zoo 


5 6+0-3 
,(,+0.09 

0.34 

or;^; 1 /on'7 
zOO.i/ zU ( 


6 2+°-® 
"•^-0.6 

19+°" 
0.34 


6 2+°° 
"•^-0.6 

Sl+°13 

0.34 

1 Ti n /01 A 
i / o.y/ zi4 


MODEL B 


kT 
Z 
Nh 

xVdof 


5 g+o-'^ 

■^■"-0.4 
22+° °^ 

16+°-^^ 

239.8/234 


0.34l°:l« 
256.1/206 


fi 5+°-" 
19+°" 

< 0.39 
180.1/208 


"-•^-0.8 
S2+°12 

< 0.97 
173.8/213 


MODEL C 


kT 
Z 

Nh 
M 

xVdof 


a £j+4.4 

6-6_i.o 

r, r,, + 0.08 

0.38l°:« 
< 2900 
238.5/233 


so+° i° 

42+°-58 
"•^"^-0.24 

< 2800 
256.0/205 


0.30t°;l^ 

< 1.29 

< 2290 
178.0/207 


1 Q+4.4 
12.3^5 g 

1-261°:^^ 

2240+^^° 
170.2/212 


MODEL D 


kT 
Z 

Nh 
M 

ANh 
X^/DOF 


r, 90 + O.O8 

"■''•J_0.08 

0.34 
1150+"°° 
< 1.1 

238.6/233 


5 4+1-3 

o,+0.08 

0.34 
< 1100 

U.C. 
255.5/205 


ll.5^«;? 
0.30l°:l^ 
0.34 

-^'^^-1340 

< 2.44 
177.9/207 


19 7+3.9 
'-6.1 
49+°-l5 

0.34 

2240+^50 

< 5.43 
170.2/212 



Notes: The best-fit parameter values and 90 per cent (Ax^ = 2.71) confidence limits from the spectral analysis of the ASCA data for 
Z-wicky 3146. Temperatures {kT) are in keV and metallicities (Z) are quoted as a fraction of the Solar value (Anders & Grevesse 1989). 
Column densities (A^h) are in units of 10^^ atom cm~^ and mass deposition rates (M) in M© yr~i . x^ values and the number of degrees 
of freedom (DOF) in the fits are given for the four spectral models discussed in Section 4.2. 



parameter in the fits. The fits to Abell 1835 and E1455+223 
in particular showed highly significant improvements with 
the introduction of this single extra fit parameter. (Note 
that absorbing material was assumed to Ue at zero redshift 
in this model). The best-fit parameter values and 90 per cent 
(Ax^ = 2.71) confidence limits from the analyses with the 
single-phase models (A and B) are summarized in Tables 
5-8. 

Although the single-phase modelling can provide a use- 
ful parameterization of the properties of the cluster gas, the 
results obtained with such a model should be interpreted 
with caution. The deprojection analyses presented in Sec- 
tion 5 show that all three of the clusters discussed in this 
paper contain large cooling fiows. The gas in the central re- 
gions of these clusters must therefore be highly multiphase 
i.e. contain a wide range of densities and temperatures at 
all radii. 

We therefore next examined the data with more sophis- 
ticated spectral models in which the spectrum of the cooling 
fiow was accounted for explicitly. The first of these models 
(Model C) consists of an RS component (to model the emis- 
sion from the ambient ICM in the region of interest) and a 
cooling-flow component (following the models of Johnstone 
e£ al. 1992) modelling the X-ray spectrum of gas cooling 
from the ambient cluster temperature, to temperatures be- 
low the X-ray waveband, at constant pressure. Note that 



Model C introduces only one extra free parameter into the 
fits relative to the single-phase Model B; the mass deposi- 
tion rate of cooling gas. The upper temperature of the cool- 
ing gas, the metallicity, and the absorbing column density 
acting on the cooling flow were tied to those of the ambient 
cluster emission modelled by the RS component. 

Fourthly, we examined a further cooling-flow model 
(D; which we expect to be the most physically-appropriate 
model) in which an intrinsic X-ray absorbing column density 
was associated with the cooling-flow. The excess absorption 
is modelled as a uniform absorbing screen in front of the 
cooling flow, at the redshift of the cluster, with the column 
density a free parameter in the fits. The column density act- 
ing on the ambient cluster emission was fixed at the Galac- 
tic value. The best-fit parameter values and confidence limits 
obtained with the multiphase, cooling-flow models (C,D) are 
also summarized in Tables 5-8. 

Finally, a fifth model in which the X-ray emission from 
the cluster was parameterized by a combination of two RS 
components was studied. However, the statistical signifi- 
cance of including the second RS component in the fits (2 
extra free parameters) is low and the results on the temper- 
atures and emission measures of the two components, which 
are poorly constrained by the data, are not presented here. 
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Table 6. Spectral Analysis of the ASCA data for Abell 1835 





Parameters 


oU 




kT 


— O.D 




z 


24+° - i2 

0.09 


MODEL A 




0.22 




X^/DOF 


416.7/369 




kT 


7 4+0.8 

6 




z 


n. ,q + 0;08 


MODEL B 




0.73+°}^ 




X^/DOF 


379.2/368 




kT 






Z 




MODEL C 


Nh 


^•"^-0.40 




M 


< 2600 




xVdof 


377.4/367 




kT 


g -r+3.5 




Z 


0.28+0:1 


MODEL D 


ATh 


0.22 




M 


200011^0 






Q 9c;+2.25 
•-'■^'-'-0.85 




xVdof 


372.9/367 



Notes: The best-fit pajrameter values and 90 per cent (Ax^ 
Abell 1835. Details as for Table 5. 

4.3 Results from the spectral analysis 

The results from the spectral analysis, presented in Tables 
5-8, provide a consistent description of the X-ray properties 
of the clusters. We find good agreement in the results from 
the different detectors. Interestingly, in a reduced sense, 
all four spectral models provide a statistically adequate de- 
scription of the ASCA spectra. Even with the single-phase 
models, however, the statistical improvement obtained by 
allowing the X-ray absorption to fit freely (i.e. the improve- 
ment obtained with Model B over Model A) is very high - 
particularly for Abell 1835 and E1455-I-223 where a simple 
-F-test (Bcvington 1969) indicates the improvement to be 
significant at >> 99.9 per cent confidence. Model B indi- 
cates excess column densities (assumed to lie at zero red- 
shift) in Abell 1835 and E1455-F223 of B.Otl'l x 10^° atom 
cm~'^and 5.9l} g x lO'^" atom cm"^, respectively. The data 
for Zwicky 3146 prefer a column density marginally less than 
the nominal Galactic value (Stark et al. 1992). 

The introduction of the cooling flow component into the 
fits with Model C results in a further reduction in x^- How- 
ever, the statistical significance of this improvement, relative 
to the model B results, is marginal, being required at > 95 
per cent confidence only with the Abell 1835 data. (Note, 
however, that the improvement obtained with model C with 
respect to model A is very high.) In general, the lowest 
values axe obtained with Model D. However, it is difficult to 
interpret the improvement in obtained with Model D with 



SI 


G2 


G3 


9.4+Ji 

0.12 


6.5+°-^ 
26+°- 12 

"■^"—0.12 


6.8+°-? 
23+°- 1? 


0.22 


0.22 


0.22 


363.5/298 


200.0/216 


199.5/261 


7 1+0.9 
7 

26+0-°^ 


6.7+?-« 
26+°' " 

'-'•^^—0.12 

< 0.82 


6.0+«-? 

23+°' " 

'-'•^^—0.10 

1 04+°''^ 


334.4/297 


199.9/215 


195.7/260 


7 9+5.0 
'•^-0.7 

0-811°:?^ 


1 2 1 +3-3 

34+0-1'' 
, Qo+0.83 


6 O+'^-i 
n 2s+° i3 

"•^•i-o.io 

1 05+10'' 
-'-■<JO_Q 52 


< 2600 


< 3500 


< 4400 


334.4/296 


198.3/214 


195.7/259 


8.011:1 

0.28«:1J 


o.32l°:l« 


"■^-0.9 
24+° l3 


0.22 


0.22 


0.22 


20001^000 


< 4300 

u.c. 


< 4300 

U.C. 


327.2/296 


198.9/214 


195.5/259 



respect to Models A-C in terms of a statistical significance 
since Model D includes fit parameters and constraints not 
present in the other models. Within the context of Model D 

the cooling flow component is required at high significance 
(Table 8). The data for AbeU 1835 and E1455-h223 also re- 
quire significant amounts of intrinsic absorption associated 
with their cooling flows (A^h ~ 3.8lJ'4 x 10^^ atom cm^^for 
AbeU 1835 and 3.8to:'^ x 10^^ atom cnr^for E1455+223) 
whereas the data for Zwicky 3146 are consistent with Galac- 
tic absorption. Note that only the SIS data have the spectral 
resolution and sensitivity at lower energies (E < 1 keV) to 
detect the presence of cooling flows in the clusters. The GIS 
data do not provide firm constraints on the emission from 
cooling gas. 

We conclude that the ASCA spectra alone are unable 
to discriminate at high significance between the multiphase 
cooling flow models (C, D) and the single-phase model (B) 
for the clusters. Although adopting Model D as intuitively 
the most reasonable description of the X-ray emission from 
the clusters leads to a strong spectral requirement for large 
cooling flows in all three systems, it is only through the 
combination of the spectral results with the results from 
the deprojection analyses discussed in Section 5, that the 
presence of massive cooling flows in these clusters is firmly 
established. 

Finally in this Section, we note the possible effects 
of uncertainties in the low-energy calibration of the SIS 
data on our results. Our own analyses of ASCA observar 



= 2.71) confidence limits from the spectral analysis of the ASCA data for 
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Parameters 


SO 


SI 


G2 


G3 




kT 
Z 


c n+0.4 
^^•^-0.3 
29+01° 


5 4+0-5 
14+0-11 


4 5+0-8 

20+0-23 
"■■^"-0.16 


4 5+0-6 

QQ+0.24 


MODEL A 




0.31 


0.31 


0.31 


0.31 




~ ,2 iT\r\T? 
X /DUi' 


178.0/170 


1 OO 1 /l OO 

138.1/138 


rt A n /oo 

94.0/83 


123.2/1U4 




kT 
Z 


4 2+0-* 
^•^-0.3 
Q ,n+0.10 


— 0.5 


4 0+1.0 

24+0-26 
0.21 


4 2+0.9 

n 41+0-25 
'-'•^^—0.20 


MODEL B 


Nh 


U.O»_(j 22 


0.89l°:^« 


< 3.0 


< 2.2 




X^/DOF 


157.7/169 


125.1/137 


92.5/82 


122.4/103 


MODEL C 


Kl 

z 


K q + 2.3 
^•"^-1.3 


5-5-1.4 

16+"" 

1 44+0-S2 


7 1 +3.1 
7.1 3 g 

34+0^25 

"■•^^-0.27 
2 5+1-2 


43+0-35 

"•^•-'_0.22 

< 3.3 




M 


< 3550 


< 3500 


< 2000 


< 2390 




X^/DOF 


155.6/168 


123.8/136 


91.6/81 


122.3/102 




kT 
Z 


c 9+2.3 
,0+0.10 


6 6+3-1 
"■"-1.9 


6.9+37 

34+0-26 


1+5.0 
0-1-1.6 
46+"-^^ 


MODEL D 


ATh 


0.31 


0.31 


0.31 


0.31 




M 




2500l^f°o 


< 2320 

u.c. 


< 2670 
U.C. 




xVdof 


154.1/168 


124.0/136 


91.4/81 


122.3/102 



Notes: The best-fit parameter values and 90 per cent (Ax^ = 2.71) confidence limits from the spectral analysis of the ASCA data for 
E1455+223. Details as for Table 5. 



tions of bright, nearby X-ray sources indicate that the cur- 
rent GSFC response matrices (released 1994 November 9) 
may slightly overestimate the low-energy response of the 
SIS instruments. This systematic effect can lead to over- 
estimates of Galactic column densities by 1 — 3 x 10^° 
atom cm~2(see also the discussion of calibration uncer- 
tainties associated with the ASCA instruments on the 
ASCA Guest Observer Facility World Wide Web pages at 
http : / /heasarc. gsfc.nasa.gov/docs/asca/caljprobs.html). 
Fixing the Galactic column densities in our analyses with 
spectral Model D at values 1.7 x 10^" atom cm~2in excess 
of the nominal Galactic values for the clusters [1.7 x 10^" 
atom cm~^being the best-fit systematic excess column den- 
sity determined from our analysis of ASCA observations 
of the Coma cluster, which we expect to contain little or 
no intrinsic absorbing material (White et al. 1991, Allen 
& Fabian 1996), we determine best-fit intrinsic column 
densities for Zwicky 3146, Abell 1835 and E1455-F223 of 
iVn < 0.6 X 10^^ atom cm"^, A^h = 3.2tl * x 10^^ atom 
cm~2, and JVh = S.llj g x 10^^ atom cm"^, respectively. 
Hence, our conclusions on the presence of intrinsic absorb- 
ing material in these clusters are essentially unaffected by 
uncertainties in the calibration of the SIS instruments. Note 
that the temperature constraints are also little-affected by 
the calibration uncertainties, with the best-fit temperatures 
and 90 per cent confidence limits from the column density- 
adjusted fits being kT = 6.9to^, 9.ltl l and 5.2^0^ keV, 
respectively. 



4.4 Single-phcise and Multiphsise temperature 
results 

The most notable differences between the results obtained 
with the single-phase (Model B) and multiphase (Models 
C,D) models are in the temperature determinations for 
the clusters. The cooling flow models imply signiflcantly 
higher ambient cluster temperatures, and therefore larger 
integrated cluster masses. 

The presence of a cooling flow will naturally lead to 
differences between the mean emission-weighted and mass- 
weighted temperatures for a cluster. The X-ray emissivity 
of the cooler, denser material in the cooling flow will be sig- 
niflcantly higher than that of the surrounding hotter gas. 
The mean emission-weighted temperature will therefore be 
biased to temperatures below the mass-weighted value for 
the system. The effects of emission-weighting become par- 
ticularly important for observations made with the compar- 
atively low spectral resolution and limited (0.1 — 2.4 keV) 
bandpass of the ROSAT PSPC. In Table 9 we present the 
results from a spectral analysis of the PSPC data for the 
central 6 arcmin radius regions of Zwicky 3146 and Abell 
1835. The single-phase models again provide a statistically 
adequate description of the spectra (xt ^ 10) but the mea- 
sured emission weighted temperatures are only 3.2+J-^ keV 
and 3.8lQ;g keV respectively, much less than 6.6+J;7 keV and 
9.51J-7 keV determined from the multiphase analysis of the 
ASCA spectra (Model D). The importance of distinguishing 
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Table 8. All instruments combined together 





Parameters 


Model A 


Model B 


Model C 


Model D 


Zwicky 3146 


kT 
Z 

Nh 
M 
X^/DOF 


5 74+0-19 
26+°-°^ 
0.34 

864.4/871 


e n'7+0.33 
26+° °* 

17+°- 12 

858.6/870 


6 9+^'^ 
"■-"-o.g 

27+ 05 
IJ-^' _0.05 

54+°-34 

1070+1270 

852.3/866 


c c+1-1 
r, 97+0.05 

< 1.02 

1330^^2?" 
853.6/866 


Abell 1835 


kT 
Z 

Nh 
M 

xVdof 


0.22 
976.4/958 


7 no+0.34 
(.UO_Q 33 

26+° °° 

72+0-11 
"•'^-0.10 

909.2/957 


^■1-1.6 
30+° 06 

-I -19+0.24 
1-1^-0.34 

898.7/953 


9.511:? 

2090l?gg 

891.8/953 


E1455+223 


kT 
Z 

Nh 
M 

xVdof 


c (11 +0.26 
^■"-^-0.26 
n. 90 + 0.06 

0.31 
541.8/501 


4 29+°-2S 

25+0-0'i' 
90+"-" 

503.5/500 


5 n+2-6 

26+°-°'^ 

1 48+°-''-''' 
1890+1300 

498.8/496 


5 4+1-9 

27+0-07 

q O + 2-0 

2030+So 
498.2/496 



Notes: The best-fit parameter values and 90 per cent (Ax^^ = 2.71) confidence limits from the spectral analyses with data from all 4 
detectors combined. Temperatures (kT) are in kcV and metallicities (Z) are quoted as a fraction of the solar Value (Anders & Grcvcssc 
1989). Column densities (A^h) are in units of lO^i atom cm~2. fcT^ z, and A^h are linked in the fits. However, the mass deposition rates 
for each detector (M; quoted in Mq yr~i ) are included as independent fit parameters, due to variations in source extraction area and 
uncertainties in the flux calibration of the instruments. The M -value quoted for models C and D is for the SO detector. 



Table 9. PSPC spectra for Zwicky 3146 and Abell 1835 



Table 10. Deprojection analyses of the clusters 



Parameter 


Zwicky 3146 


Abell 1835 




^cool 


^cool 


M 


kT 
Z 

Nh 

xVdof 


o 9+1-4 
•^•^-0.7 
n 59+1-68 
"•^^-0.48 

23+0-05 
26.3/22 


q + 1.6 

•^•°-0.9 
o+°-22 

n 1O+0.03 

U.1S_„ Q2 

17.0/20 


Zwicky 3146 

Abell 1835 

E1455+22:-! 


1 01+°-°° 

-1 00+0.31 
l-"^''-0.31 

1-1' (i.i.", 


231^50 
231+2° 


1OC5+408 
1106lf25 

--.,,,-:i(,2 

' (i4 











Notes: The best-fit parameter values and 90 per cent (Ax^^ = 
2.71) confidence limits from the spectral analysis of the PSPC 
data. 

between single-phase and multiphase models in the analy- 
sis of X-ray data for clusters is discussed in more detail in 
Section 7.1. 



5 DEPROJECTION ANALYSIS 
5.1 General results 

Wo have carried out a deprojection analysis of the ROSAT 
images using an updated version of the code of Fabian et al. 
(1981). Using assumptions of spherical symmetry and hy- 
drostatic equilibrium in the ICM, the deprojection technique 



Notes: A summary of the results from the deprojection analyses 
of the ROSAT HRI data. Cooling times (< cooi ) are mean values 
for the central (8 arcscc) bin and arc in units of 10^ yr. Cool- 
ing radii (rcooi), the radii at which the cooling time exceeds the 
Hubble time (1.3 X lO"' yr), are in kpc. Integrated mass deposi- 
tion rates within the cooling radii (M ) are in units of Mr yr^i . 
Errors on the cooling times arc the 10 and 90 percentile values 
from 100 Monte Carlo simulations. The upper and lower confi- 
dence limits on the cooling radii are the points where the 10 and 
90 percentiles exceed, and become less than, the Hubble time, 
respectively. Errors on the mass deposition rates are the 90 and 
10 percentile values at the upper and lower limits for the cooling 
radius. Galactic column densities as listed in Table 1 are assumed. 
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Figure 7. A summary of the results from the deprojection analysis of the HRI data for Zwicky 3146. Prom left to right, top to bottom, we 
plot; (a) surface brightness, (b) pressure, (c) integrated luminosity, (d) temperature, (e) electron density, (f) cooling time, (g) integrated 

gas and gravitational mass and (h) integrated mass deposition rate. Data points arc mean values and Icr errors (in each radial bin) from 
100 Monte Carlo simulations, except for (d), (f) and (h) where the median and 10 and 90 percentile values have been plotted. 



can be used to study the properties of the intracluster gas 

{e.g. density, pressure, tonipcraturo, cooling rate) as a func- 
tion of radius. The deprojection method requires that either 
the total mass profile (which defines the pressure profile) or 
the gas temperature profile be specified. Following ASCA 
observations of nearby cooling flow clusters (Fabian et al. 
1996), and the results from the combined X-ray and gravi- 
tational lensing study of the cooling-flow cluster PKS0745- 
191 (Allen et al. 1996), we assume that the mass- weighted 
temperature profiles in the clusters remain constant, at the 



temperatures determined from the fits with the multiphase 

spectral models to the combined detector data sets (Table 
8). Spectral model D is intuitively the preferred model, but 
using the full temperature range allowed by models C and D 
probably provides a more realistic estimate of the true un- 
certainty on the mass-weighted temperatures in the highly 
complex, multiphase environments of the cooling flows. Col- 
umn densities were are flxed at the Galactic values from 
Stark et al. (1992), but see also Section 5.3. 

The clusters discussed in this paper are remarkably 
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Figure 8. A summajry of the results from the deprojection analysis of the HRI data for Abell 1835. Details as for Fig. 7 



similar in their X-ray properties to PKS0745-191 and the 
analogy to that system is a reasonable one. It should be 
noted that although the deprojection method of Fabian et al. 
(1981) is essentially a single-phase technique, it produces re- 
sults in good agreement with the more detailed multi-phase 
treatment of Thomas, Fabian & Nulsen (1987) and, due to 
its simple applicability at large radii, is better-suited to the 
present project. 

The azimuthally averaged X-ray surface brightness 
profiles of the clusters determined from the HRI data 
(background-subtracted and corrected for telescope vi- 
gnetting) and the results from the deprojection analyses are 



summarized in Figs. 7—9. The primary results on the cooling 
flows in the clusters; the central cooling times, the cooling 
radii and the integrated mass deposition rates within the 
cooling radii are listed in Table 10. The results on the cool- 
ing flow in Zwicky 3146 are in good agreement with those 
reported by Edge et al. (1994) from an earlier analysis of 
the HRI data. The mass deposition from the cooling flows is 
distributed throughout the inner ~ 200 kpc of the clusters 
with M J5 r. As noted in Section 1, such distributed mass 
deposition profiles requires that the central ICM is inho- 
mogeneous (Nulsen 1986; Thomas, Fabian & Nulsen 1987; 
Fabian 1994). Note also that the mass deposition profiles 
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Figure 9. A summary of the results from the deprojection analysis of the HRI data for E1455+225. Details as for Fig. 7 



shown in Figs. 7 — 9(h) flatten at radii < 200 kpc, where the 
cooling time is < 10^° yr. Thus accounting for the look-back 
time to the clusters (~ 4 x 10® yr) does not significantly 
alter the integrated mass deposition rates. 



the X-ray gas mass/total mass ratios for the central 500 
kpc of the clusters (where the total mass is assumed to be 
described by the best-fit parameters listed in Table 11). 



5.2 Parameterization of the cluster masses 

In Table 11 we summarize the mass distributions required to 
produce the flat temperature proflles described in Section 5. 1 
We have parameterized the mass distributions as isothermal 
spheres (Equation 4-125 of Binney & Tremaine 1987) with 
core radii, rc, and velocity dispersions, a. In Fig. 11 we plot 



5.3 The correction for intrinsic absorption 

As discussed in Section 5.1, for the purposes of deprojection 
wc have assumed that the column densities to the clusters 
are given by the Galactic values of Stark e£ al. (1992). How- 
ever, the spectral analysis of Section 4 clearly shows that 
both Abell 1835 and E1455-I-223 exhibit significant excess 
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Figure 4. (Upper Panel) The SO spectrum for Zwicky 3146 with 
the best-fit solution for Model B overlaid. The data have been 
binned-up by a factor 5 along the energy axis for display purposes. 
(Lower Panel) Residuals to the fit. The positive residuals at i? ~ 
0.55 keV are due to small systematic uncertainties in response 
matrix around the oxygen edge in the detector. 



Figure 5. (Upper Panel) The SO spectrum for Abell 1835 with 
the best-fit solution for Model B overlaid. (Lower Panel) Residuals 
to the fit. Details as in Fig. 4. 



Table 11. Cluster mass distributions 



kT 



Zwicky 3146 

Abell 1835 
E1455-I-223 



6.6 
9.5 
5.4 



-1-3.0 
0.7 
-1-4.9 
2.0 
-1-2.2 
1.1 



(li:?) 
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10001?°° 
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Notes: A summary of the temperature constraints (in keV) and 
the velocity dispersions (in kms~^ ) and core radii (in kpc) of the 
isothermal mass distributions (Binncy & Trcmainc 1989) required 
to produce the flat temperature profiles shown in Figs. 7 — 9(d). 
Errors on the temperatures give the the full range allowed by 
spectral models C and D. (The tighter constraints determined 
with model D alone are given in parentheses.) Errors on the ve- 
locity dispersions show the range of values required to match the 
temperature results. 
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Figure 6. (Upper Panel) The SO spectrum for E1455-I-223 with 
the best-fit solution for Model B overlaid. (Lower Panel) Residuals 
to the fit. Details as in Fig. 4. 
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Figure 10. The ratio of tlie gas mass to the total mass as a 
function of radius. 



absorption. This can have a significant effect on the mass 
deposition rates determined from the dcprojcction analyses. 

To correct the deprojection results for intrinsic absorp- 
tion, we adopt spectral model D eis the most reasonable de- 
scription for the clusters. The intrinsic column densities in 
Abcll 1835 and E1455+223 (iVn = 3.8 x 10^^ atom cm"^; 
Tabic 8) are assumed to be in uuiform screens in front of 
the cooling flows. The effects of absorption on the observed 
HRI count rates from the cooling flows have been calcu- 
lated from XSPEC simulations, using the ROSAT HRI re- 
sponse matrix issued by GSFC. For both Abcll 1835 and 
E1455-I-223, the intrinsic absorption acts to reduce the count 
rates from the cooling flows by a factor of two (in detail, 
factors 2.07 and 2.04, respectively). The true mass deposi- 
tion rates from the cooling flows can therefore be assumed 
to be a factor 2 larger than the values listed in Table 10. 
[For Zwicky 3146, the maximum allowed intrinsic column 
density of 10^^ atom cm~^implies a maximum correction 
factor to the mass deposition rate of 1.28.] The absorption- 
corrected mass deposition rates for the clusters arc summa- 
rized in Table 12. The agreement between these values and 
the spectrally-determined mass deposition rates (Model D) 
is excellent. 



6 OPTICAL PROPERTIES OF THE CLUSTERS 
6.1 The galaxy populations 

The galaxy populations in Abell 1835 and E14554-223 have 
been studied using the Palomar B and I images of the clus- 
ters. The SExtractor software of Bertin (1995) was used to 
identify galaxy candidates from their total (Kron) I mag- 
nitudes and {B — I) colours (measured in 3.0 arcsec di- 
ameter apertures from seeing-matched images). For Abell 



1835, 283 objects were selected with 15.0 < I < 21.0 and 
2.5 < (B - /) < 3.2. The CCG (which is bluer than al- 
lowed by this range) was also included, making 284 galaxy 
identifications in total. For E1455-f223, 188 objects with 
15.0 < / < 21.0 and 1.8 < (B - /) < 2.7 were identified. 
Absolute V magnitudes were calculated from the observed I 
magnitude, and applying the appropriate K-corrections. 

In Figs. 11 (a), (b) we plot the (projected) galaxy pro- 
files and the ratio of the galaxy mass to the total mass for 
Abell 1835 and E1455+223. A value for (M/L)v of 10 has 
been cissumed. Note the large central peaks in the ratio pro- 
files [Fig. 12(b)] due to the luminous central galaxies, and 
the near-flatness of the ratios at large radii (r > 0.5 Mpc) 
indicating that the galaxies follow an approximately isother- 
mal distribution. 



6.2 Spectra of the central cluster galaxies 

Optical spectra for the CCGs of Zwicky 3146 and Abell 1835 
are presented and discussed by Allen (1995). These CCGs 
are two of the most optically line-luminous central galaxies 
known, with luminosities in HaA6563 emission alone of ~ 
10'*^ ergs~^ . In Fig. 13 we present the spectrum of the 
CCG in E1455-I-223. The data were obtained with the Faint 
Object Spectrograph on the Isaac Newton Telescope (INT) 
in June 1991. A slit width of 1.5 arcsec was used, providing 
a resolution of ~ 16 A in the first spectral order. 

The CCG of E1455+223 exhibits strong, narrow emis- 
sion lines (FWHM ~ 500 kms"'^) and an enhanced blue 
continuum with respect to a typical elliptical galaxy spec- 
trum. Both features are characteristic properties of CCGs in 
large cooling fiows (Johnstone, Fabian & Nulsen 1987; Heck- 
man et al. 1992; McNamara & O'Connell 1989; Crawford & 
Fabian 1993; Allen 1995). The observed flux in HaA6563 



of 1.9 ± 0.1 X 10 org cm ^s ^implies a slit luminosity 
of 6.1 ± 0.3 X lO"*^ ergs"^ . However, the H/3A4861 emis- 
sion line is only marginally detected and, after correcting 
for absorption by the underlying stellar continuum (a fac- 
tor ~ 2 correction), we derive an HaA6563/H/3A4861 flux 
ratio of > 4.0. This implies significant intrinsic reddening 
at the source [E{B — V) > 0.3] and an intrinsic HaA6563 
(slit) luminosity of > 10**^ ergs"'^ . [The reddening laws of 
Seaton (1979) and Howarth (1983) have been used.] Don- 
ahue, Stocke and Gioia (1992) present results from a narrow- 
band Hq-|-[NII] imaging study of E1455-I-223. These authors 
show that the line emission from the CCG is highly extended 
and that the total line flux exceeds the slit flux reported here 
by a factor ~ 7 — 8. The Donahue et al. (1992) results, to- 
gether with the spectral results reported hero, thus imply an 
HaA6563 luminosity for the CCG of E1455+223 of > 7x 10"*^ 
ergs"'^ , comparable to that of Zwicky 3146 and Abell 1835. 

The CCGs of Zwicky 3146 and Abell 1835 also exhibit 
significant intrinsic reddening. In Table 12 we summarize the 
results on intrinsic X-ray absorption and optical reddening 
for the clusters. These results, together with the optical/X- 
ray/UV results discussed by Allen et al. (1995) for a larger 
sample of cooling flows at intermediate redshifts {z ~ 0.15), 
indicate an interesting tendency for clusters with large col- 
umn densities of intrinsic X-ray absorbing material to ex- 
hibit significant intrinsic reddening. The similarity of the 
column densities inferred from the X-ray, optical and UV 
data, across a variety of aperture sizes, suggests a dust-to- 
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Figure 12. Optical spectrum of the CCG of E1455+223 obtained with the Faint Object Spectrograph on the INT in June 1991. 



Table 12. Reddening and Excess absorption 







Spectral M 


Corrected deproj. M 


E{B - V) 


.'Vh.opt 


Zwicky 3146 
Abell 1835 
E1455+223 


< 1.0 

o + 1.6 
3-O-0.4 
o O+2.0 
3-8-0.8 


20901700 


1355±SI 

22911^3? 
149ll?30 


22+°- IS 

49+° !^ 
"■^^-0.15 

> 0.3 


1 O+0.9 
9 S+1-0 

^ 1-7 



Notes: Columns 2 and 3 list the intrinsic X-ray column densities (without account for systematic uncertainties in the SIS calibration; 
Section 4.3) and mass deposition rates determined from the spectral analysis of the combined instrument data sets (Table 8). Column 4 
lists the mass deposition rates, determined by deprojection, corrected a posteriori for the effects of X-ray absorption (Section 5.3). For 

Zwicky 3146, zero intrinsic absorption has been assumed. Applying the maximum allowed correction factor for this cluster (1.28) gives a 
mass deposition rate of 1734+20(3 ■ E(B — V) estimates in column 5 are determined from the Ha6563/H/34861 line ratios in the CCGs. 
The data for Zwicky 3146 and Abell 1835 are from Allen (1995). A'^h.opt values are the column densities of X-ray absorbing material 
implied by the E{B — V) estimates, following the relation of Bohlin, Savage & Drake (1978). 



gas ratio in these galaxies similar to that in our own Galaxy. 
The results also suggest that much of the dust may be asso- 
ciated with, or entrained within, the X-ray absorbing gas. 



7 DISCUSSION 

We have discussed, in detail, the X-ray properties of Zwicky 
3146, Abell 1835 and E1455+223. Wo have shown that all 
three of these clusters contain exceptionally large cooling 
flows. Zwicky 3146 and Abell 1835 are amongst the most 
X-ray luminous clusters known (Table 2). The cooling flows 



in these systems account for ~ 15 — 20 per cent of the total 
intrinsic luminosity in the 2—10 keV band, and as much 
as ~ 40 per cent in the 0.1 - 2.4 keV ROSAT band. With 
E1455+223, which is a factor 2 — 3 less luminous than the 
other clusters, the cooling flow accounts for ~ 35 per cent 
of the luminosity in 2-10 keV band and ~ 60 per cent of the 
emission between 0.1 and 2.4 keV. 

Both AbeU 1835 and E1455-I-223 exhibit significant in- 
trinsic absorption in their ASCA spectra. The need for ex- 
cess absorption is found in both the single-phase and multi- 
phase (cooling flow) spectral analyses and cannot reasonably 
be attributed to uncertainties in the Galactic column den- 
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sities (Stark et al. 1992). The most plausible interpretation 
of the excess absorption is that it is due to material asso- 
ciated with the cooling flows (spectral Model D). The mass 
of absorbing gas implied by an intrinsic column density of 
~ 3.8 X 10^^ atom cm , distributed in a uniform (circular) 
screen across the central 220 kpc (radius ~ rcooi) of the clus- 
ters, is ~ 4.6 X 10^^ Mq. Such a mass could plausibly be ac- 
cumulated by the cooling flows in Abell 1835 and E1455-I-223 
in ~ 2 - 3 X 10^ yr. [See also White et al. (1991) and Allen 
et al. (1993).] Note also the excellent agreement in the mass 
deposition rates for the cooling flows determined with the 
spectral and deprojection methods (Section 5.3) under this 
assumption for the distribution of absorbing gas. With spec- 
tral models B and C, wherein the excess absorption is as- 
sumed to cover the whole cluster, the mass of absorbing gas 
is implausibly high. The intrinsic column densities inferred 
for Abell 1835 and E1455-I-223 are similar to those observed 
in nearby cooling flows (White et al. 1991; Allen et al. 1993; 
Fabian et al. 1994; Fabian et al. 1996). Note also that the 
metallicities of Z ~ 0.25 — O.SOZq measured for these clus- 
ters are similar to those observed in nearby systems, and 
imply that the the bulk of the enrichment of the ICM in 
these clusters occurred before redshifts of ~ 0.3. 




0.5 1 

Radius (Mpc) 



Figure 11. (a) The projected galaxy mass as a function of radius 
in Abell 1835 and E1455-I-223. (b) The ratio of the galaxy mass 
to the total mass (derived using the temperature constraints from 
spectral Model D). 



7.1 Cooling flows and multiphase models 

One of the most important results from this paper is the 
marked difference in the temperatures of the clusters deter- 
mined from the single-phase and multiphase spectral mod- 
els. These results are summarized in Table 13. The single- 
phase temperatures consistently (and significantly) under- 
estimate the multiphase results. For Zwicky 3146 and Abell 
1835, the ASCA single-phase results underestimate the mul- 
tiphase temperatures by ~ 10 and 25 per cent, respectively. 
With the ROSAT data, the discrepancy in much more se- 
vere, with the PSPC values underestimating the ASCA mul- 
tiphase results by ~ 3.4 keV (50 per cent) for Zwicky 3146, 
and ~ 6 keV (60 per cent) for Abell 1835. Note, however, 
that in all cases the reduced values indicate statistically 
acceptable fits. 

The multiphase kT results (Models C,D) should approx- 
imate the true mass-weighted temperatures in the clusters 
(Thomas et al. 1987; Waxman & Miralda-Escude 1995; Allen 
et al. 1996). The single-phase results, however, are simply 
emission/detector- weighted average values for the integrated 
cluster emission. Since the X-ray emissivity of cluster gas 
rises with increasing density (decreasing temperature), the 
presence of a large cooling flow naturally leads to a decrease 
in the emission-weighted temperature of a cluster. The ef- 
fects on the emission-weighted kT are most dramatic in the 
0.1 — 2.4 keV band of the PSPC, where the emission from 
cooler gas phases dominates the detected flux. However, the 
(comparatively) poor spectral resolution and limited band- 
pass of the PSPC, moan that the single-phase models can 
still provide a statistically adequate description of the data. 
The PSPC data are unable to discern the need for multi- 
temperature components (although the imaging data clearly 
require them). 

These results imply that caution should be applied in 
the interpretation of temperatures determined with simple, 
single-phase models and, in particular, those determined 
from ROSAT data. With ASCA data the single-phase re- 
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Table 13. Comparison of multiphase and single-phase results 



PSPC (SP) ASCA (SP) ASCA (MP) 

Cluster kT xl (y) kT xl {^) kT xl {^) 

Zwicky 3146 1.20 (22) e.llg 3 0.987 (870) 6.6lg ^ 0.986 (866) 

Abell 1835 3.8t Jig 0.85 (20) 7.Qto% 0.950 (957) 0.936 (953) 



Notes: A comparison of the temperature results obtained with the Single-phase (SP) and MultiPhase (MP) models. SP results from 
ASCA are for spectral Model B. MP results are for spectral Model D. 



suits should be more reliable, although significant discrep- 
ancies can still arise (as in the case of Abell 1835). 

The presence of a range of density and temperature 
phases is clearly established by the data for cooling flow 
clusters. However, it should not assumed that the absence 
of a cooling flow implies that a single-phase modelling of 
the ICM is appropriate. The existence of cooling flows with 
distributed mass deposition requires significant inhomogene- 
ity (a density /temperature spread of ~ a factor 2) in the 
ambient cluster gas before the cooling flow forms (Nulsen 
1986; Thomas, Fabian & Nulsen 1987). The best data for 
clusters are consistent with such a range of inhomogeneity 
(Allen et al. 1992b). The absence of cooling flows in some 
nearby, luminous clusters such as the Coma cluster is usu- 
ally attributed to merger events having disrupted the cluster 
cores and having re-heated and redistributed the cooling gas 
throughout the cluster. In such circumstances it seems un- 
likely that the merger will completely homogenize the gas 
and, therefore, that a single-phase model will provide an 
exact measure of the mass-weighted cluster temperature. 

7.2 A comparison with lensing masses 

In Fig. 13 we show the Palomar U band image of the central 
2.5 arcmin^ of Abell 1835. The CCG is the bright source in 

the centre of the field. The image shows a number of dis- 
torted features (arcs, arclets and image pairs) attributable 
to gravitational lensing by the cluster. In particular, we ob- 
serve a bright, elongated arc ('A' in Fig. 13) at a radius of 
30.3 arcsec from the centre of the CCG, along a PA of 133 
degree. The arc has a length of ~ 16 arcsec. is extended 
along a PA of 221 degree, and exhibits reflection symmetry 
about the point 14''01'"03.7s, 02°52'21". 

For a simple, circular mass distribution the projected 
mass within the tangential critical radius, ret, is given by 

Mp.oj(rct) = ^ (jT-ir' ) 

y -^clus-^arc— clus y 

where ret can be approximated by the arc radius (rare = 
150 kpc) and -Ddus, i'arc and -Darc-cius are respectively the 
angular diameter distances from the observer to the clus- 
ter, the observer to the lensed object, and the cluster to the 
lensed object. In Fig. 14 we show the mass within arc 'A' 
as a function of the rcdshift of the arc, calculated with the 
above formula. Also shown are the X-ray constraints on the 
projected mass within this radius [1.0 x 10^* Mq< Mproj < 
2.1 X 10^* M© from the multiphase analysis using spectral 



Model C, and 1.1 x 10^"* Mq< Mp,oj < 1.6 x 10^^ Mq using 
spectral Model D]. Combining the X-ray and lensing mass 
results we axe able constrain the redshift of arc 'A' to be 
> 0.7. It is also important to note that if the single-phase 
(Model B) X-ray temperature results for Abell 1835 were 
(wrongly) used in place of the multiphase values, no con- 
sistent solution for the X-ray and lensing masses would be 
possible. 

The intrinsic ellipticity of the lensing potential may 
lead to a slight (< 20 per cent) overestimate of the lens- 
ing mass determined with the circular mass model (Bartle- 
mann 1995). However, the effects of ellipticity also lead to 
a slight overestimate of the X-ray mass within this aperture 
and, to first order, the conclusions on the redshift of the arc 
should not be dramatically affected. [The lensing properties 
of Abell 1835, and those of a larger sample of X-ray lumi- 
nous clusters observed with the Palomar 5m telescope, are 
discussed further by Edge et al. (1996). ] 

Small et al. (1995) report results from a study of 
(weakly) gravitationally distorted images in the field of 
E1455-I-223, from which they derive a projected mass within 
450 kpc of the cluster centre of ~ 3.6 x lO" Mq. This 
mass exceeds the X-ray determination of the projected mass 
within this radius, 1.6 x 10^* Mq (for an isothermal mass 
distribution corresponding to a temperature of 5.4 keV), by 
a factor ~ 2. [Note that the determination of the X-ray mass 
assumes that the cluster remains isothermal and extends to 
3Mpc. However, extrapolating the mass profile to 5 Mpc 
increases the projected mass within 450 kpc of the cluster 
centre by only 2 per cent.] Using a potential consistent 
with the upper-limit to the X-ray temperature of 7.6 keV 
{i.e. a = 900 kms"^ , rc = 40 kpc) we still determine a 
projected mass within 450 kpc of only 2.5 x 10^* Mq. 

The lensing result on the cluster mass for E1455-I-223 
appears high. E1455-I-223 is a regular, relaxed cluster with 
a large cooling fiow, and a 2 — 10 keV X-ray luminosity of 
1.3 X 10**^ erg s~^ . The ASCA constraints on the X-ray tem- 
perature (4.3 — 7.6 keV) are consistent with results for other 
nearby, cooling-flow clusters of similar X-ray luminosity {e.g. 
AbeU 1795; Edge et al. 1990, Fabian et al. 1996), which lends 
support to the X-ray mass determination. The velocity dis- 
persion of cr = 660 - 900 kms"^ (Table 11) implied by the 
X-ray data is also in good agreement with optical obser- 
vations (cr ~ 700 kms~^ ; Mason et al. 1981, Smail et al. 
1995). 

A cluster of exceptional X-ray luminosity and tem- 
perature is required to provide a projected mass within 
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Figure 13. The Hale 5m U band image of Abell 1835. Arc 'A' is indicated to the South East of the CCG. 



450 kpc consistent with the Smail et al. (1995) result for 
E1455+223. Abell 1835, discussed in this paper, provides a 
meiss of 2.5 — 5.3 x 10^'' M©. Similarly, the exceptionally 
X-ray luminous cooling-flow cluster PKS0745-191 (Lx = 
2.8 X 10*^ ergs"^ ), for which Allen et al. (1996) present a 
self-consistent determination of the mass distribution from 
X-ray and gravitational lensing data, provides a projected 
mass within 450 kpc of only ~ 3.7 x 10^"* M©. Given the 
X-ray luminosity of E1455+223 (which is a factor 2-3 less 
than PKS0745-191 or Abell 1835), the X-ray mass measure- 
ment for the cluster seems reasonable and the lensing mass 
high. The result of Smail et al. (1995) may imply an un- 
usual redshift distribution for the weakly distorted sources, 
or a projected mass distribution that deviates significantly 
from the simple isothermal mgiss model used (perhaps due to 
some line of sight mass enhancement from material external 
to the X-ray luminous part of the cluster) . 



7.3 Optical and X-ray properties 

CCGs in cooling flows frequently exhibit characteristic low- 
ionization, optical emission-line spectra {e.g. Johnstone, 
Fabian & Nulsen 1987; Heckman et al. 1989; Crawford & 
Fabian 1992; Allen et al. 1995). The optical (HqA6563) line 
luminosity correlates with the excess UV/blue continuum lu- 
minosity (the excess with respect to the UV/blue emission 



expected from a normal gE/cD galaxy; Johnstone, Fabian 
& Nulsen 1987; McNamara & O'ConneU 1989; Allen et al. 
1992a; Crawford et al. 1995; Allen 1995). Typically, both the 
emission lines and the excess UV/blue continua are extended 
across the central 10 — 20 kpc of the clusters. 

The clusters discussed in this paper contain three of 
the largest cooling flows known. They also host three of 
the most optically line- luminous (and UV/blue luminous) 
CCGs. Only the CCG of the massive cooling flow cluster 
PKS0745-191 exhibits a comparable optical line luminos- 
ity (Allen et al. 1996). Although no simple correlation be- 
tween M, tcooi and the optical line luminosity (and therefore 
UV/blue continuum luminosity) exists, the data presented 
here confirm a tendency for the most optically-line-luminous 
CCGs to be found in the largest cooling flows {e.g. Allen et 
al. 1995). 

The UV/blue continua in Abell 1835 and Zwicky 3146 

appear dominated by emission from hot, massive stars 
(Allen 1995). These young stellar populations may also pro- 
vide the bulk of the ionizing continuum emission responsible 
for the observed optical emission lines. The exceptionally 
high mass deposition rates from the cooling flows in the 
clusters can naturally provide the large reservoirs of cooled 
material necessary to fuel the observed (very high) star for- 
mation rates (Allen 1995). The star formation may also ac- 
count for some of the dust in the CCGs. 
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Figure 14. A comparison of the X-ray and lensing mass estimates 
for arc 'A'. Tlic bold curve sliows tlic mass within the arc deter- 
mined from tlic standard lensing formula (for a splierical mass 
distribution) as a function of the rcdsliift of the arc. The solid 
horizontal lines show the X-ray constraints on the mass within 
this radius determined with spectral Model D (7.8 < kT < 10.8 
keV). The dashed lines show the constraints for spectral Model C 
(7.5 < kT < 14.4 keV). The lensing and multiphase X-ray results 
together imply Zarc > 1-6 for Model D, and Zarc > 0.7 for Model 
C. The lower dotted line shows the mass within arc 'A' implied 
by the single-phase spectral results (Mproj 9.1 X 10-*-^ MqIot 
kT = 7.0 keV). The single-phase X-ray results are inconsistent 
with the lensing data. 

The excellent alignments between the (optical) CCG 
and (X-ray) cluster isophotes are consistent with the results 
for other massive cooling flows at intermediate {z ~ 0.15) 
redshifts (Allen et al. 1995). These results again reveal the 
unique and intimate Unk between CCGs and their host clus- 
ters. 



8 CONCLUSIONS 

Wo have presented detailed results on the X-ray proper- 
ties of Zwicky 3146, Abell 1835 and E1455+223, three of 
the most distant, X-ray luminous clusters known. We have 
shown that these clusters contain the three largest cooling 
flows known, with mass deposition rates of ~ 1400, 2300 
and 1500 M© yr~^ , respectively. We have presented mass 
models for the clusters and have highlighted the need for 
multiphase analyses to consistently explain the spectral and 
imaging X-ray data for these systems. The inappropriate use 
of single-phase models leads to significant underestimates of 
the cluster temperatures and masses. For Aboil 1835 it was 
shown that a mass distribution that can consistently explain 
both the X-ray and gravitational lensing data for the cluster 
can only be formed when multiphase X-ray spectral models 



are used. We have discussed the relationship between intrin- 
sic X-ray absorption and optical reddening in the clusters. 
These results suggest that the X-ray absorbing material fre- 
quently observed in the X-ray spectra of cooling fiows is 
dusty. 
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